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Normal muon decay, n + — > e + u e u IM , is studied as a tool to discriminate between 
the Dirac and Majorana types of neutrinos and to survey the structure of the weak 
interaction. It is assumed that massive neutrinos mix with one another and that the 



interaction Hamiltonian consists of the V — A and V + A charged currents. A new 
set of parameters used in place of the Michel parameters is proposed for the positron 
distribution. Explicit forms of these new parameters are obtained by assuming that 
the masses are less than 10 eV for light neutrinos and sufficiently large for heavy 
Majorana neutrinos, which are not emitted in the muon decay. It is shown that a 
possible method to discriminate between the Dirac and Majorana cases is to use a 
characterization given by the x 2 fitting of their spectra. It is also confirmed that the 
theoretical predictions in the Majorana neutrino case are almost the same as those 
obtained from the standard model. Indeed, their differences cannot be distinguished 
within the present experimental precision. 
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§1. Introduction 



The structure of the leptonic charged weak interaction provides us with an important 
source of information that may lead to a unified theory beyond the standard model. Normal 
muon decay is a pure leptonic process accessible to precise measurements of this structure 
with high statistics, because it is free from the complications of the strong interaction and 
hadronic structure. 

Experimental data have been analyzed by employing the helicity-preserving four fermion 
weak interaction with (S ± P), (V ± A) and T forms, ^ because this arrangement allows one 
to make direct contact with specific models. The Michel parameters have been used to obtain 
some information concerning the structure of the weak interaction under the assumption that 
the neutrinos are massless and the lepton number is conserved. Recent experimental data 
have exhibited smaller deviations from the predictions based on the standard model. 2 ^ 

The neutrino emitted in the annihilation of negatively charged leptons has been regarded 
as a particle, while the neutrino created together with the negatively charged leptons has 
been assigned to an anti-particle. In this assignment, in which there is a distinction between 
the neutrino particle and its anti-particle, the neutrino is of the Dirac type, and the lepton 
number is conserved in the weak interaction. In the other case, in which there is no such 
distinction, the neutrino is of the Majorana type. In the Majorana case, the lepton number 
is not conserved. 

The dominant interaction responsible for muon decay has a V — A structure, and the 
standard model is constructed on this footing. In this model, the left-handed neutrino field 
v L is assigned to a member of a doublet of the SU(2) L x U(l) group, and no right-handed 
field v' R is present. The neutrino is assumed to be massless [m v = 0), and the charged 
current weak interaction takes place via the exchange of the left-handed weak gauge boson 
Wl- The neutrino and anti- neutrino have the definite helicities h = —1/2 and h = +1/2, 
respectively, and we cannot distinguish between the Dirac and Majorana neutrino within 
the standard model. 

Now, however, it has been established through the discovery of neutrino oscillation that 
neutrinos have finite masses and mix with one another. 5 ^ Thus, it is now known that the 
neutrino cannot be in a definite helicity state. For this reason, it is possible to discriminate 
between neutrinos of the Dirac and Majorana types. It is an important and fundamental 
question to determine whether the neutrino is of the Dirac or Majorana type. 

There are two other unsolved problems regarding leptons. One is to determine why the 
observed mass differences among the three neutrinos are so small in comparison with charged 
leptons and quarks. The other is to understand why the left-handed V — A interaction is 
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favored over the right-handed V + A interaction which has not been detected definitively. In 
the framework of gauge theory, it seems natural that the V — A interaction is favored, as a 
result of the spontaneous breakdown of left-right symmetry, which is believed to be satisfied 
at sufficiently high energy. 

One appealing way to solve these problems simultaneously is to use the idea of the seesaw 
mechanism, through which the right-handed neutrino field v' R is introduced. 6 ^ Let us explain 
the scenario of this mechanism for the one generation case for simplicity. In this scenario, the 
terms related to the neutrino mass can have the Majorana-type mass term (u' r ) c Mrv' r , in 
addition to the Dirac-type mass term, v' r Mdv^. The Majorana neutrino field is defined after 
diagonalizing the mass matrix, 7 ) (see Eq. (IA-3j) in Appendix A). The left-handed Majorana 
neutrino can have a small mass, of the order of M^/Mr, while the right-handed Majorana 
neutrino can have a large mass, of the order of M R , provided that the condition <C M R 
is satisfied. This situation leads to various extensions of the standard model, such as the 
SU{2) L x SU{2) R x U(l) model. If this scenario is indeed valid, then we can conclude that 
there exists a right-handed gauge boson Wr and that the muon decay receives contributions 
from interactions whose structure is somewhat different from that of the standard model. 

It is not yet known whether the neutrino is of the Dirac or Majorana type and, fur- 
ther, what structure the weak interaction has beyond the standard model. For this rea- 
son, it is necessary to construct a method that provides some information concerning these 
points. Neutrinoless double beta decay which violates lepton number conservation is the 
only presently known possible way to directly determine the type of the neutrino. However, 
this decay process requires very high resolution experimentally, because of the long half-lives, 
due to the tiny neutrino mass and/or the small contribution from the V + A current. Thus, 
this method is yet to provide a decisive conclusion. 

Muon decay takes place irrespective of the type of neutrino involved and precise data 
have been accumulated and analyzed to investigate the structure of the weak interaction 
by assuming the neutrino to be massless and of the Dirac type. It is also meaningful, as a 
complementary study, to survey the possibility of whether muon decay can be used as a tool 
to determine the type of the neutrino. With this in mind, we note that there is a difference 
between the spectra of emitted positrons in the Dirac and Majorana neutrino cases. 

The aim of this paper is to propose a new parameterization of muon decay that is suitable 
for analyzing the type of neutrinos and the structure of the weak interaction. We adopt 
a Hamiltonian consisting of both V — A and V + A currents, which is inspired by the 
SU{2)l x SU(2)r x £7(1) gauge model, and present a method to analyze the implications of 
the experimental data. 

In 321 we summarize the general framework of our study and discuss kinematical effects 
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on the emitted due to finite neutrino mass (m„ 7^ 0) . There, our assumptions and 
the approximation adopted in our analysis are discussed. In £0 the e ± energy spectrum is 
surveyed in detail. We propose some parameterizations to discriminate between the types of 
neutrinos and discuss their experimental feasibility. The polarization of the e ± is discussed 
in 21 A summary and conclusion are given in £0 In Appendix A, to make the paper 
self-contained, features of the lepton mixing matrix for a model with left- and right-handed 
neutrinos are summarized. Also there, the details of the coupling constants for the weak 
interaction Hamiltonian based on the SU(2)l x SU(2)r x U(l) gauge model are presented 
for convenience. In Appendix B, definitions of various coefficients are listed, and their 
explicit forms under certain conditions are given for the Dirac and Majorana neutrino cases 
separately. 

§2. General framework 

We assume the following form of the effective weak interaction Hamiltonian for the /x ± 

decay: 7 ) 

Kw(x) = ^{ji La J« L + \jl R J1 R + ml Ra ^ L + kjIlJZb) + H.c. , (2-1) 

where Gf is the Fermi coupling constant. The left-handed and right-handed charged weak 
leptonic currents, jn and jm, are defined as 

3£La(x) = £ e (x)j a (l - 7 5 )^lO) and jm a {x) = £i{x)^ a {l + ^)v' m {x). (2-2) 

Here £%, v tL and v tR are the weak eigenstates of the charged lepton, left-handed neutrinos 
and right-handed neutrinos, with flavors I = e and fi. The interaction in Eq. (j2-lj) is a general 
form of the four fermion, derivative-free Lorentz-invariant interaction, which consists of the 
V — A and V + A currents. 

The weak eigenstates of the charged leptons (£1) and neutrinos [yg L and Vm) are defined, 
respectively, as the superpositions of the mass eigenstate charged leptons Eg and neutrinos 
Nj with mass rrij. Then, the charged currents are expressed as 

2n 

jeLa(x) = y^J5/(z)7 a (l - i5)UtjNj(x), 

3 2n (2-3) 

jm a (x) = ^ Ei(x)j a (l + 7 5 )%iVj (a) , 
3=1 

for the case of the n generations. 7 - ) Here Ugj and Vjy are, respectively, the left-handed and 
right-handed lepton mixing matrices. 
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The weak interaction given in Eq. (|2T|) is naturally expected from the gauge models that 
contain the left-handed and right-handed weak gauge bosons, Wl and Wr. In these models, 
the appearance of the coupling constant A is due to Wr, while terms with 77 and k come 
from the possible mixing between Wl and Wr. As a typical example with right-handed 
interactions, we consider the SU(2)l x SU(2)r x £7(1) model in Appendix 1X1 If this model 
is assumed, then the coupling constants k and 77 in Eq. (|2-lj) can be taken as identical, as 
shown in Eq. (|AT9J) . However, they are treated as independent constants in this paper in 
order to allow comparison with the more general case without a restriction from the gauge 
theory (see, e.g., Ref. 1)). The structures and magnitudes of the lepton mixing matrices U 
and V are also briefly summarized in Appendix [XJ We do not take account of the mirror 
lepton currents and Higgs boson exchange, for simplicity. 

Now we study the normal muon decay /i + — ► e^vJT^ (or \i~ — > e~T^u^). In the framework 
of the effective weak interaction Hamiltonian given in Eq. ([2-ip . ^ decay takes place as 



/i ± - e ± + Nj + N k , (24) 

where N k represents an antineutrino for the Dirac neutrino case, but it should be understood 
as Nk for the Majorana neutrino case. 

If the radiative corrections are not included, 8 '' the differential decay rate for polarized e ± 
in the rest frame of polarized ^ is expressed as 1 )' 9 -* 

Loose ^-4AD ix ,e)ii + p eix ,e).a (2-5) 

where 

E ml + ml 

x = — and W = — ^ = 52.8 MeV. (2-6) 

W 2m M v ' 

Here m^ and m e are the muon and electron masses, respectively, and E is the energy of e ± . 
The angle 9 is the direction of the emitted e ± with respect to the muon polarization vector 
at the instant of the /j ± decay. By taking account of the finite neutrino mass (m u 7^ 0), 
the allowed range of x is limited kinematically as 

x ^ x ^ x max = ( 1 - rj fc ), (2-7) 

where 



x = = 9.65 • 10- 3 and r% = V J = 3.63 ■ 1(T 14 . (2-8) 

Here mj and mfc are masses of the neutrinos emitted in the muon decay and should satisfy 
the relation 



rrij + m k < (m M - m e ). (2-9) 
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In Eq. ()2-8j) . the neutrino masses have been taken as rrij = m k = 10 eV, in order to obtain a 
rough idea of the magnitude of r 2 k . 

The constant A in Eq. (j2-5J) is introduced to simplify the expression for the differential 
decay rate by using the arbitrariness of its normalization. It is referred to as a normalization 
factor in this paper. It is shown in ^21 that there are various possibilities for the choice of A, 
when experimental data are analyzed, although these choices differ only by rearrangements 
of the terms in the theoretical expression. This constant A plays the role of the leading term 
in the overall normalization related to the muon lifetime. 

In the differential decay rate given in Eq. ()2-5|) . D(x, 8) is the energy spectrum part, 
expressed as 

D(x, 6) = [N(x) ±P tl cos6P(x)\, (2-10) 

where P M = |P M |, and the functions N(x) and P(x) are, respectively, the isotropic and 
anisotropic parts of the e ± energy spectrum. Their details are discussed in JJJ The plus 
(minus) sign in Eq. ()2-10|) corresponds to n + (/x - ) decay. The vector P e (x, 0) in Eq. ()2-5|) 
is a polarization vector of e ± , and £ is the directional vector of the measurement of the 
spin polarization. We discuss P e {x, 6) in 

The isotropic and anisotropic parts of the energy spectrum consist of various terms, 
each of which has a different x dependence taking a complicated form and includes some 
combination of the lepton mixing matrices. Therefore, in order to extract useful information, 
it is desirable, as the first step, to investigate the characteristic features of their x dependences 
and to simplify them. 

As an example, let us consider only two terms of N(x), 

x 2 ) + (b + - a+) (2x - x 2 - xl)] , (2-11) 

b + = (b + X 2 b), (2-12) 
^ = ^jfc FjC^jk {Uej^lU^kl 2 , (2-13) 
b = t 3k F) j2 G 3k \V e] \ 2 \V ilk \ 2 . (2-14) 

Here, the primed sum represents the sum taken over neutrinos whose emission in the /x ± 
decay is allowed by the restriction given in Eq. (|2-9J) . 

Two kinematical factors, Fj k and Gj^, emerge from the phase space integral of the 
emitted neutrinos in both the Dirac and Majorana neutrino cases. They represent the 
additional kinematical effect due to the finite mass of the neutrino {m v ^ 0), and their 
explicit forms are given in $ 12.11 [see Eqs. (|2-15j) and ([2 -16)) ]. The treatment of the lepton 
mixing matrices, Uij and V^j, is considered in ^2.21 



N(x) = 




where 

a + = (a + A 2 a), 

a = ^jk F jk 2 \ U ej\ 2 \U^k 

a = T! jk F% 2 \V e3 \ 2 \V, k \ 
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2.1. Kinematical factors due to the non-vanishing neutrino mass ■■■ Condition (A) 

First, it is worthwhile to note that all terms in the spectrum D(x, 9) and the polarization 

— * — - 1/2 /i ■• \ 

P e (x, 9) in Eq. (|2-5|) are proportional to F-' k , for example, as shown in Eq. (|2-13|) . These 
terms include some combinations of the following kinematical factors:*** 

„2 



jk 



r 



j k 



G jk - 

where fijk is defined by 



1 — x 

2 



r 



j k 



X 



1 - 



1 + 2 



1 — X 

2 



r 



.1 k 



1 



X 



1 -Afijk) 
(!- fJ-jk) 



+ 6 



'jk 



X 



H'j ki 



(2-15) 
(2-16) 



mj mk 



(rrij +m k y 



(2-17) 



If m v = 0, both Fjk and Gjk are unity for the entire range of x, as seen from Eqs. ()2-15j) 
and ()2-l(ij) . However, if m v ^ 0, they display the following different types of behavior near 



X rl 



jk 



and G 



J k 



6/i 



jk 



in the limit x 



x T1 



1 — r 



(2-18) 



This implies that the spectrum and polarization of e tend to zero suddenly near x r 



Therefore, in principle, the shape of the spectrum near x T 
the massive and massless neutrino cases. 



exhibits different behavior for 



Fortunately, the x dependences of Fjk and Gjk appear significantly only in a very tiny 



range near x max , say, for x > (1 — 10 8 r 



jk) 



(1 — 10 6 ) if the required numerical accuracy of 
the experiment is of order 10~ 6 . Thus, practically, there seems to be no problem in treating 
*) Strictly speaking, some terms proportional to ra v include another kinematical factor accompanied by 

1/2 

Fj k (see Ref. 10)). However, we do not consider these terms in this paper, because their contributions are 
negligible. 

**' Here we state the reason why Fjk suddenly becomes zero near a; max for m v ^ 0, in spite of the fact that 
Fjk = 1 for m v = 0. It is convenient to introduce the momentum transfer squared, A 2 = (qj + qk) 2 , where 
qj is the 4-dimensional momentum of Nj. Because the e ± energy E is given by E — (m 2 ^ + m 2 — A 2 )/2m lll 
the maximum energy is realized when A 2 takes the minimum value A 2 = (m,j + mk) 2 , as shown in Eq. I|2-7|l . 
This minimum is realized under the following two conditions: (i) Both of the neutrinos are emitted in the 
direction opposite to e in order to satisfy the momentum conservation, namely, q~j + qi + p e = 0. (ii) 
Each neutrino has a definite momentum, say qj = p e mj/(rrij + rrik)- Therefore, there is no freedom to 
assign arbitrary neutrino momentum. In other words, the density in the phase space is zero; that is, we 
have Fjk — 0. A similar situation exists in the case that the mass of one neutrino is zero. Contrastingly, if 
vrij = m,k = 0, we have to choose A 2 = 0. This situation is allowed for various combinations of two neutrino 
momenta, since only the total momentum of the neutrinos is fixed under the condition (i). Thus, there is no 
special restriction in the phase space compared to the general case for an arbitrary e energy; that is, we 
have Fjk = 1. 
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Fjk and Gj k in D(x, 9) and P e {x, 9) as independent of x. More specifically, if we assume that 
the emitted neutrinos have small masses, i.e., at most of the order of 10 eV, the following 
approximations yield very good accuracy: 

F jk = 1 and G jk = 1. (2-19) 

Hereafter, this approximation is referred to as Condition (A). 

Under Condition (A), the second term of N(x) in Eq. (|2-llj) yields no contribution. 
Hence, we need not consider this term, except for the negligible tiny range near x max . We 
omit such terms in this paper without any loss of practical accuracy. 

2.2. Contributions from the lepton mixing matrices U and V ■ ■ ■ Condition (B) 

The sum over the square (or product) of the lepton mixing matrix elements appears in 
the spectrum D(x, 9) and polarization P e {x, 9) appearing in Eq. ()2-5|) . for example as seen 
from Eqs. (j2-lHj) and (j2-14)l . If Condition (A) is accepted, we are able to take the sum over 
neutrino indices under some assumptions. 

In the Dirac neutrino case, it is assumed that all neutrinos can be emitted in the fi^ 
decay. Then we have the following properties from the unitarity conditions of U and V: 

E j \U ej \ 2 = 'Z j \Vt j \ 2 = l. (2-20) 

By contrast, in the Majorana neutrino case, we assume the existence of heavy Majorana 
neutrinos, which are not emitted in the // ± decay. Then, there is a different situation, in 
which we have 

Sjl^f = 1 - ui 2 and ^-\V ej \ 2 = !T e 2 , (2-21) 

where the primed sums are taken over only the light neutrinos. Here ui 2 and W 2 are the 
representatives of small deviations from unitarity due to these heavy Majorana neutrinos. [ 
The details are given in Eqs. (jA-6|) - (jA-8|) . (jA-15|) and (jA-16|) .] In what follows, the explicit 
forms of these matrices U and V are not needed. 

In addition, in the Majorana neutrino case, the following products of U and V appear: 

= Ej U ej Vpj and w^Th = V ek U^ k . (2-22) 

An example in which w efl appears is mentioned in ^2 .HI The quantities w eil and w^Th are 
also small, as shown in Appendix |X] and Eq. (|4-14|) . The assumptions Eqs. ()2-20|) - ()2-22|) 
are referred to as Condition (B). 
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2.3. Terms characteristic of the Majorana-type neutrino 

Finally, let us explain why the Majorana-type neutrino offers different information from 
the Dirac-type neutrino. We emphasize that some contributions that are specific to the 
Majorana neutrino case are added to the decay rate. Among these, one appears even in the 
(V — A) interaction alone if m v ^ 0. The others appear if the (V + A) interaction exists in 
addition to (V — A), even in the case m u = 0. We study two examples. 

As the first example, let us consider the mode in which both the (e + , Nj) and (fi + , Nk) 
vertices are of the (V — A) type. This assignment is referred to as the ordinary mode (A). If 
m v ^ 0, Nj has a small component of helicity h = +1/2, whose magnitude is proportional 
to ruj/uj, where ujj is the energy of Nj. Therefore, in the Majorana neutrino case, in which 
this Nj cannot be distinguished from the Nj associated with there is a new cross mode 
(B) for which Nj (= Nj) comes from the annihilation of fi + and Nk (= Nk) is emitted with 
e + . The interference between the ordinary mode (A) and the cross mode (B) appears in the 
decay probability. It becomes proportional to the product of the neutrino masses, mjUik- 
Of course, there is no such possibility in the Dirac neutrino case, where Nj ^ Nj. 

Now, if the (V + A) interaction is added to (V — A), the situation changes greatly. As a 
second example, let us consider the mode for which the (e + , Nj) vertex is of the (V — A) type, 
but the (/i + , N),) vertex is of the (V + A) type. In this case, the helicity of Nk is h = —1/2, 
the same as that of Nj, even in the limit m u —>■ 0. We call this case the mode (C). In the 
Majorana neutrino case, it is possible to have another cross emitting mode (D) with the (e + , 
Nk) and (fi + , Nj) vertices, in addition to the mode (C). There arises interference between 
the (C) and (D) modes in the decay probability. The leading term of this interference is not 
proportional to m u , but it includes the coupling constant k 2 with lepton mixing matrices 
\^jUejVnj\ 2 . This results from the equality of the helicities, mentioned above. The mode 
(D) does not exist in the Dirac neutrino case. [For details, see Fig. (2) and Table I of Ref. 9)]. 

Many terms proportional to m u exist in the decay rate in both the Dirac and Majorana 
neutrino cases, because of the small component of the helicity that is proportional to m„. 
However, these terms are negligibly small and are not taken into account in this paper. The 
complete decay formulae including these terms are given in Refs. 10) and 11). 

§3. Energy spectrum of 

The isotropic part N(x) and anisotropic part P(x) of the e energy spectrum appearing 
in Eq. (|2T0|) are, respectively, expressed as follows under Conditions (A) and (B) given in 
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N(x) = [a + (3x - 2x 2 - x 2 ) + 12 ( k +c + e m k +m )x{l-x) 

+Qe m \d r XQ (1 - x)} , (3-1) 

P(*) = (£) \i* 2 -4 + ^ - rl) 

+ 12(k_ c + e m k_ m )(l-x)]. (3-2) 

Here, the decay formulae for the Dirac and Majorana neutrinos are obtained by setting 
e m = and e m = 1, respectively. The original forms of N(x) and P(x) before adopting 
Conditions (A) and (B) are presented in Eqs. ()B-2jl and ()B-3j) of Appendix B. 

The first terms of N(x) and P(x) represent their x dependences obtained from the stan- 
dard model, where A = a± = 1, and all other coefficients (k± c , k± m and d r ) are zero. The 
quantity Tq in P(x) is defined as follows: 

Before listing the details of these coefficients, we note here that the suffix c of k± indicates 
contributions that are the same for the Dirac and Majorana neutrino cases, whereas the suffix 
m indicates coefficients associated with the Majorana neutrino only. 

In the Dirac neutrino case, these coefficients are expressed as follows: 



a±= (1± A 2 ) and k± c = {k 2 ± rj 2 ) . 



(3-4) 



Of course, in this case there are no contributions from k± m and d r . Note that the final 
expressions for N(x) and P(x) in the Dirac neutrino case are the same as those obtained by 
assuming m v = for one generation if Conditions (A) and (B) in §2] are both satisfied. 
In the Majorana neutrino case, these coefficients have the following complicated forms: 

a±= [(l-u- e 2 )(l-u; 2 )±X 2 - e 2 v^ 2 ], (3-5) 



k± c = U « ^i-«r)v ±v % (i-vO > ( 3 - 6 ) 



k± m = Qj [k 2 I w eti | 2 ±rf | w eiih | 2 ] 



(3-7) 



d r = Qj Re(w eM * w eiMh ). (3-8) 

Here, We 2 , Ue 2 , and w efl h are small quantities in general, as mentioned with regard to 
Eqs. ()2-21|) and ()2-22|) . It should be noted that k± c has the same order of magnitude as 
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k± m , in contrast to the Dirac neutrino case. It is worthwhile noting that the final N(x) and 
P(x) in the Majorana neutrino case exhibit small deviations from those of standard model 
independently of the magnitudes of the coupling constants A, k and r\. 

Since 1985, experimental data have been analyzed using expressions based on the well- 
known Michel parameters. Expressing N(x) and P(x) in terms of the Michel parameters 
using our notation, we have 

r 2 1 

N(x) = 6 x(l - x) + - p M (4x 2 — 3x — Xq) + r] M x (1 - x) , (3-9) 

9 

(l-x) + ~5 M (4x-3-r 2 )]. (3-10) 



P(x) = 2f M \ h 2 -x 2 



These expressions are obtained, for example, from Ref. 1) or Eqs. (31) and (32) of Ref. 12).*) 
Hereafter, these forms are referred to as the Michel parameterization. They are presented 
for the Dirac neutrino case with m v = 0. In the standard model, these parameters take 
definite values: pu = 5m = 0.75, £m = 1, and rj M = 0. In our model, there is no r] M term, 
even in the massive Dirac neutrino case, if we ignore terms proportional to m„. The reason 
why the rjM ^ term appears in the Michel parameterization of Ref. 1) is that it comes 
from the interference between the (V ± A) and (S ± P) (or T) interactions. 

At first glance, the relation between N(x) in Eq. (|3-lj) and that in Eq. ()3-9j) may not be 
clear. It is shown in ^3.11 that the Michel parameterization given in Eq. ()3-9|) is a special 
case of our expression appearing in Eq. ()3-l|) in which some appropriate constant value is 
chosen for the normalization factor A. A similar situation is demonstrated for P(x) in £13.21 

3.1. Isotropic part of the spectrum: N(x) 

First, let us introduce the following quantity for the normalization factor A: 

A n£ = (a+ + 2nk +c + 2£e m k +m ) > 0, (3-11) 

where n and I are some integers. Next, we verify that Eq. (|3-9J) is a special case of Eq. ()3-l|) 
with the choice of either A = A 10 ot An. 

We first consider the case with A = Aw, in which terms characteristic of the Majorana 
neutrino case appear explicitly. Then, the following expression is derived: 



N(x) 



(lb) {^ a + + 2k +c)^-2x 2 -xl) 



+A; +C [12 x{l-x)- 2(3x - 2x 2 - x 2 )] 



*' Our definitions and those in Ref. f) are related as N(x) = 6Fjg(x) and P(x) = 6Fas{%)- The 
amplitudes g]^ defined in Ref. 1) correspond to our coupling constants as follows: g\ L = 1, g^ R = A, 
9lr = K ' 9rl ~ 7 h with all other amplitudes set to zero in the present paper. 
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+e m [12fc +m x(l-x) + 6\d r x (l-x)]} (3-12) 
(3a; - 2x 2 - x 2 ) + 2 p c (3s - 4a; 2 + x 2 ) 

+ 12e m p m x(l- x) +6e m r] m x (l- x) . (3-13) 

Here, the parameters are defined as follows: 

Pc = 0*) >0 , Pm . (^) > o and Ifm = g£) . (3,4) 

The fact that p c and p m are positive is clear from Eqs. ()3-4|) - (|3-7j) and (j3-ll|) . Now, it is 
easy to confirm that the Michel parameterization, Eq. ()3-9|) . can be obtained from Eq. ()3-13|) 
by introducing the relation 

2p c = (i-^PMj , (3-15) 

because, for this confirmation, it is not necessary to take account of the terms p m and r) m 
appearing only in the Majorana neutrino case. Note that the relation A = 1 is required in 
the Michel parameterization [see Eq. (44) of Ref. 12)]. 

Next, if the case with A = An is chosen, we again obtain the form given in Eq. ()3-9|) if 
we replace p c in Eqs. ()3-13|) and ()3-15j) by (p c + e m p m ), because the third term in Eq. ()3-13|) 
is absorbed into the second term through this replacement, and the fourth term plays the 
role of the tjm parameter. Of course, the denominator v4 10 of Eq. (|3-14j) should be replaced 
by A 1V 

Thus, it is clear that there is no advantage of introducing the Michel parameter at 
present, although it was very useful in determining the type of the weak interaction. In 
fact, in the Michel parameterization defined in Ref. 1), the form {pM — 3/4), like Eq. ()3-15j) 
[and (5m — 3/4), like Eq. ()3-20|) ] is used. Therefore, instead of measuring the deviation from 
Pm = 0.75, it seems desirable to directly determine p c and p m themselves, which indicate 
the deviations from the standard model. Furthermore, it is worthwhile to note that we are 
interested in estimating the x dependence of the difference between the experimental data 
and the prediction of the standard model, namely, (3a; — 2x 2 — Xq) appearing in Eq. (jHHJ), 
through use of the x 2 -fitt m g i n the experimental analysis. 

Recently, the TWIST group reported a precise experimental result for Pm^ 

p M = 0.75080 ± 0.00032(stat) ± 0.00097(syst) ± 0.00023, (3-16) 

where the third error comes from the ambiguity in t]m appearing in Eq. ()3-9|) . This ambiguity 
is due to the fact that various values for t)m have been used within the uncertainty of the 
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accepted average value r)M = (—7 ± 13) • lO -3 . 1 ** Assuming the Dirac-type neutrino, this 
group obtained the result |tan£| < 0.030 by combining Eqs. (|3T5|) . ()3T4|) and ()3-4|) with 
Eq. (jA-19|) . where ( is the Wl-Wr mixing angle defined in Eqs. (jA-13|) and (jA-14|) . 

Although the TWIST group analyzed their data using Eq. (j3-9|) . their results can be 
interpreted in our analysis as follows. If neutrinos are of the Dirac type, there is no difference 
between the choices A = Ai and A = An, because in this case p m = r\ m = 0. For the 
Majorana neutrino case, we choose A = An. Then, the restriction on (p c + e m p m ) can be 
obtained from their result for pu by using the relation in Eq. (j3-15|) . and some information 
regarding T] m can be extracted from their interpretation of the t]m term. 

It is useful to note that we have p^ < 0.75 in our model with the choice A = An for 
either neutrino type. This follows from the relation (p c + e m p m ) > 0, as seen in Eqs. ()3-14|) 
and ()3-15|) . But the mean value in Eq. ()3-16|) is pu = 0.75080. Accordingly, this implies 
(p c + £ m Pm) < 0, although p M < 0.75 is satisfied within experimental uncertainty. *> 

Thus, we cannot distinguish from this experimental result whether the neutrino is of the 
Dirac or Majorana type. Although the r] m parameter in Eq. ()3-13|) is characteristic of the 
Majorana neutrino, this parameter cannot be used for this distinction, because not only is 
it multiplied by the small coefficient Xq, but also rj m itself takes a very small value in our 
model. 

However, there is a method that might make it possible to distinguish between the two 
neutrino types. This method makes use of the different x dependences of the coefficients for 
p c and p m , as seen from Eq. (|3-13|) . For example, suppose we analyze experimental data by 
using Eq. ()3-13|) with p m ^ and obtain some x 2 value, say, x m f° r the Majorana neutrino 
case. Then, suppose we repeat a similar analysis using Eq. ()3-13|) with p m = i] m = and 
thereby determine \\ for the Dirac neutrino case. Indeed, such a Xd was already determined 
by the TWIST group. 2 - ) In any case, if Xm ^ s m uch smaller than x% we can conclude that 
there is a higher probability that neutrinos are of the Majorana type. 

As an aside, it is worthwhile noting that p c and p m appear symmetrically. For example, 
we have the same coefficient for k± c and k± m in Eq. (I3TJ1 . which corresponds to the choice 
A = Aqq. Thus, the choice A = A nn is not of interest to us, because it does not discriminate 
between the Dirac and Majorana neutrino cases. On the other hand, if we choose A = A i, 
we have an expression similar to Eq. ()3-13|) for A = Aiq. That is to say, the roles of p c and 
p m are exchanged. Of course, theoretically, A\q and A i are different, but it is not easy 
to distinguish them experimentally from the total decay rate, because their deviations from 
unity seem to be small. 

*) In order to account for the old result, pm = 0.7518 ± 0.0026, larger values of \gf L \ and \g% L \ have 
been obtained in Ref. 1). 



13 



We now discuss an important property useful for choosing A when we try to compare 
Xm with Xd- ft i s desirable to make the difference between Xm an d x\ as l ar g e as possible. 
Because this difference depends on the different x dependences of terms including the p c 
and p m parameters, this situation can be realized by choosing A n £ with n ^ £, as seen in 
the derivation of Eq. ()3T2|) . Let us fix £ = for simplicity. Then, the x dependence of the 
term including p m is 12x(l — x), as shown in Eq. ()3T3j) . while the x dependence of the term 
including p c is listed for various n in Table I. Therefore, we have the freedom to minimize 
the x 2 value by choosing A n g according to the pattern of data distribution. Despite this 
fact, hereafter, we use the choice A = A 10 in this and subsequent sections to simplify our 
description. 

Table I. The x dependence of the term including p c for the cases of various A n Q. 



n 




Term including p c 





A) 


12x(l - x) 


1 


A 10 


2x(3 - Ax) + 2x 2 


2 


A 20 


-Ax 2 + Ax 2 


3 


A 30 


— 6x + 6xq 


4 


^4 


-4x(3 -x) + 8x 2 



3.2. Anisotropic part of the spectrum: P(x) 

Next let us examine the features of P(x) and introduce the following quantity for the 
common factor: 

B ni = (a- + 2nk„ c + 2£e m k„ m ). (3-17) 

We now show that Eq. ()3-10p is a special case of Eq. ()3-2|) with the choice of either B 30 
or B 33 . First, in the case of B 3Q , we obtain the expression 

P(x) = ^x 2 -x 2 {(-l +2x- r 2 Q ) + 6 5 C (3 - Ax + r 2 ) + 12 e m 5 m (1 - x)}, 

(3T8) 

where the parameters are defined as follows: 

It is easy to obtain Eq. (|3T(J|) from Eq. (|3T8j) by introducing the relations 

6 6 C = (l - ~ 5 M \ and £ = Cm, (3-20) 
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because we do not take account of the term 5 m appearing only in the Majorana neutrino 
case. 

Next, in the case with -B33, we obtain the very same form as Eq. ([310)) . with S c in 
Eqs. (13-18)) and ()3-2(J|) replaced by (5 C + e m 5 m ). Of course, the constant B 30 in Eq. ()3-19j) 
should be replaced by B 33 . 

It is worthwhile noting that 5 C and 5 m also appear symmetrically. Therefore, the situation 
is quite similar to that with the choice of the normalization factor A = A n £. It is possible 
to obtain different x dependences of the terms including 5 C and 5 m by choosing a common 
constant B n £ with n 7^ I. By again fixing i = for simplicity, the x dependence of the term 
including 6 C is tabulated for various values of n in Table II. However, this property is not so 
effective in this case of P(x), in contrast to the case of N(x), because 5 C and S m themselves 
become zero or very small, as discussed in the next paragraph. 

Table II. The x dependence of the term including 5 C for the cases of various B n Q. 



n 


B n 


Term including S c 





Boo 


12(1 -x) 


1 


Bio 


2(7 -Sx + rl) 


2 


B20 


4(4 - 5x + rl) 


3 


B30 


6(3-4x + rg) 


4 


B40 


4(5 - 7x + 2rl) 



If the SU(2) L x SU{2) R x U{1) model is assumed, then the relation 



K — Tj (3-21) 

is obtained from Eq. (jAT9|) . Therefore, we can conclude that in the Dirac neutrino case, we 
have 

S c = 0, (3-22) 

as seen from Eqs. f)3- 19j) and ()3-4|) . On the other hand, in the Majorana neutrino case, the 
parameters (5 C and 5 m ) include some very small contributions coming from the differences 
between small values of the lepton mixing matrices, as seen from Eqs. ([3 -6)) and (|3-7|) . Thus 
in this case, although strictly speaking they are non-zero, for practical purpose we can 
approximate them as zero: 

5 C ~Q and 5 m ~ 0. (3-23) 
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From these relations, we can conclude within the SU(2)l x SU(2)r x U(l) model that 
we have the following expression: 

P{x) = ^x 2 - xg(-l + 2x - r 2 ). (3-24) 

If we set £ = 1, this P(x) itself is that obtained from the standard model. Although this form 
is common to the Dirac and Majorana neutrino cases, it should be noted that the definition 
of £ differs in the two cases, as seen from Eqs. (j3-19|) . (|3-4j) - (|3-6|h namely, 

(1 — A 2 ) 

£ = t — — for the Dirac neutrino case, (3-25) 

(l + y + 2 rj z ) 

£ ~ 1 for the Majorana neutrino case. (3-26) 

Here, the deviation from £ = 1 for the Majorana neutrino case is not expressed explicitly, 
because it cannot be measured within the present experimental precision. This definition of 
£ is independent of the choice of B n £ in this SU{2) L x SU(2) R x U{\) model. 
Note that the Michel parameter 5m should be 

5 M = 0.75, (3-27) 

within the present experimental precision, as seen from Eq. (|3-2(Jj) . because 5 C ~ in this 
model. The values 5m = 0.75 and £ = 1 seem to be consistent with the recent experimental 
results obtained by the TWIST group; 3 ) 

5 M = 0.74964 ± 0.00066 ± 0.00112, (3-28) 
0.9960 < P M £<£< 1.0040. (3-29) 

Under the assumption that neutrinos are of the Dirac type, this group estimated A c < 
(80.425/420) 2 = 3.7 • 10~ 2 from Eqs. (j3~2U|l and (j3~2K|) . where A c is approximately equal to 
the ratio of the mass of Wl to the mass of Wr squared [see Eq. (IA-18|) ]. 

3.3. Summary for the spectrum 

Let us summarize our expression for the spectrum. Because it consists of many terms, we 
ignore some small terms, like xo and Tq, in order to see the essential features. Furthermore, 
in order to see its characteristic features, we assume the SU(2)l x SU{2)r x U(l) model; in 
other words, we ignore the 5 C and 5 m parameters in Eq. (|3-18|) . Then, we have the expression 

D(x, 6)=x [{3 - 2x) + 2 p c (3 - 4x) + I2e m p m {l - x) ± cos6 {-1 + 2x)] , 

(3-30) 
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where the normalization factor A = A\ $ has been used. If another x dependence of the term 
including the p c (or p m ) parameter is chosen according to Table I, all A±o in this section 
should be replaced by the corresponding A n £. 

Finally, we mention the theoretical expression for the experiment which determined the 
following quantity: 

'P(x) 



oj = P., lim 

M X— 1 



N{x) 



(3-31) 



The experimental result 

uj > 0.99682 (3-32) 

was reported by Jodidio et al. 13 ^ 

Here, it should be noted that the real allowed range of x is limited by x max = (1 — r| fc ), 
as shown in Eq. (J2-7j) . Also, both N(x) and P(x) become zero at x max , as mentioned 
in Eq. (|2-18|) . However, since these restrictions are only effective in a very tiny range, 
it is understood that this definition of u is a theoretical result obtained by taking the 
extrapolation from the allowed range of x. In addition, because the radiative correction is 
known to be larger for x > 0.9, as shown in Ref. 8), it is assumed that the experimental 
results are adjusted by taking this radiative correction into account. 

Because this u is defined by taking the ratio of two parts of the spectrum, it is 
independent of the choice of the normalization factor A n g. In other words, its theoretical 
expression is obtained from Eqs. (|3-1|) and (|3-2|) as follows: 

u = P li — = Pf, , (3-33) 

a + pm 

where the coefficients a± are defined in Eqs. ()3-4j) and (|3-5|) . and the last expression is 
obtained from the Michel parameterization given in Eqs. ()3-9|) and ()3-10j) . For simplicity, 
the contribution from the electron mass is not included here by setting x^ = Tq = 0. 

For convenience, we now present explicit forms of uo in our model. In the Dirac neutrino 
case, it is 

<" = P m^> (3-34) 
while in the Majorana neutrino case, it becomes 

(1 W/) (1 - Iv 2 ) + A 2 -» 2 " ^' {6 6b} 

Here, the last approximation is good in practice, because ug 1 and 2 seem to be very small, 
as shown in Eq. (|2-21j) . In other words, no deviation from the standard model can be 
expected again in the Majorana neutrino case, within the present experimental precision. 
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§4. Polarization of e 1 * 1 

We now define three components of the spin polarization of e ± , P e (x, 9) in Eq. (2-5). 
Its longitudinal component along the momentum direction (p e ) is expressed as Pl{x,9). In 
order to separate its transverse components, we choose the decay plane defined by this p e and 
the muon polarization vector The components of the transverse polarization within 

and perpendicular to this decay plane are, respectively, expressed as Pti{x, 9) and Pt2{x, 9). 
Mathematically, these three components are expressed as follows: ^ 

P e (x, 9) = P L {x, 9)p e + P T1 (x, g) (P« x P,) x p, + Pt2(X; Q) ■PeXP, (4-1) 

\{Pe X Pfj,) Xp e \ \p e X P^l 

It is convenient to separate the x-dependent parts of these components from the emission 
angle of e ± , namely, cos 9 = (p e ■ P^j. Therefore, we introduce the quantities 

p / m _ ±Q(x)+PyCos9S(x) 

D(x,9) ' 1 j 

where the denominator D(x, 9) is given in Eq. ()3-30|) . if the SU{2) L x SU{2) R x U{\) model 
is assumed. 

The explicit expressions of these Q(x), S(x), R(x) and T(x) are presented in terms of 
the parameters defined in §3. Strictly speaking, if A = A n £ is chosen in the case of N(x), 
Ai appearing in this section should be replaced by A n £. 

4.1. Longitudinal polarization: Q(x) and S(x) 

The isotropic part, Q(x), and anisotropic part, S(x), of the longitudinal polarization are, 
respectively, expressed as follows under Conditions (A) and (B) given in ^21 



Q(x) =^x*-xt[(3-2x- rl) - 6 4(1 - r 2 ) + 12 e m 5 m (1 - x)] , (4-5) 
S(x) = [(-x + 2x 2 -x 2 ) + 2p c {7x-8x 2 + x 2 ) 

+12 e m p m x(l - x) -2e m 7] m x (1 - x)] . (4-6) 

The parameters (^, 5 C and 5 m ) in Q(x) are defined in Eq. ()3-19j) for the case of P(x). In other 

words, the common factor B30 has been used.*-* In addition, the parameters (p c , p m and i] m ) 

*' This function Q(x) is equal to 6Fjp(x) of Ref. 1), where one new parameter, £' M , was introduced, 
because the common factor Bio was used instead of our -B30, and the interference term between the (S±P) 
and T interactions was included in £' M in order to obtain a compact expression. Therefore, Fip(x) has a 
different x dependence. 
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in S(x) are defined in Eq. ()3-14j) for the case of N{x).*> Of course, the corresponding results 
from the standard model are obtained by setting £ = 1 and all other parameters to zero. 
The original forms of Q(x) and S(x) before adopting Conditions (A) and (B) are presented 
in Eqs. (JB~H) and (JH-5}) of Appendix B. 

If the SU (2)l x SU {2)ji x U (1) model is assumed, Q(x) is reduced to the following simple 
expression corresponding to Eq. (|3-24|) : 



Q(x) = £y/x 2 -x 2 (3-2x-r 2 ). (4-7) 

Thus, the longitudinal polarization is expressed as follows, if some small terms, like Xq, Tq, 
5 C and 5 m , are ignored in order to see the essential features: 

-P„ cos [(1 - 2 x) - 2 p c (7 - 8.i) -l2e mPm (l- x)\}. (4-8) 

Therefore, the longitudinal polarization of e ± from an unpolarized muon (or 6 = tt/2) is 
expressed as follows in the SU{2) L x SU(2) R x U(l) model: 



P L (x, tt/2) = ±£ 



2p c (3-4x) 



for the Dirac neutrino, (4-9) 



(3-2a;) + 2p c (3-4a;)_ 
Pl{x, tt/2) = ±£ for the Majorana neutrino, (4-10) 

where £ is given in Eqs. (J3-25|) and (J3-26j) for the respective cases, and the additional terms 
due to the parameters p c and p m have been ignored for the Majorana neutrino case, because 
of their smallness in comparison with the present experimental precision. 

Burkard et al. 14 -' reported the following experimental result for e + from an unpolarized 
muon by assuming p c = for the Dirac neutrino case: 

P L (x, tt/2) = 0.998 ± 0.045. (4-11) 

The present average value given in the Particle Data Group is Pl = 1-00 ± 0.04. 1 ) 

4.2. Transverse polarization within the decay plane: R(x) 
Under Conditions (A) and (B) presented in ^21 we have 

R(x) = [-(1 - 14 p c - 12 e m p m ) x Q {l - x) - 2 e m rj m (x - x 2 )] . (4-12) 

*' This function S(x) is equal to 6Fap(x) of Ref. 1), where one new parameter, was introduced, 
because the common factor A^q was used instead of our Am, and also the interference term between the 
(S ± P) and T interactions was included in Therefore, Fap(x) has a different x dependence than our 
S(x). 



19 



The result of the standard model is obtained by setting p c = p m = r] m = 0. The original 
form before adopting Conditions (A) and (B) is presented in Eq. (jB-61) .*- 1 

In the Dirac neutrino case, we cannot expect any useful information from this measure- 
ment, because the first main term obtained from the standard model is already proportional 
to the small Xq. By contrast, the Majorana parameter r\ m is accompanied by the coefficient 
(x — Xq), which is significantly larger than the coefficient x (l — x) in N(x) [cf. Eq. (|3T3j) ]. 
However, r\ m itself is small, as discussed in the next paragraph. Therefore, it is rather difficult 
to obtain definite information regarding 7] m from this measurement. 

Let us estimate the order of magnitude of the parameter r) m , which includes the following 
combination of lepton mixing matrices: 

Xd r = X Re£; fe (UZjVekV*!/^) = X Re(w^ * w^K), (4-13) 

where Eqs. (jSUl), flUBl) and have been used - 

Since it can be naturally assumed that there are no contributions from the heavy Majo- 
rana neutrinos, we are able to express as follows by omitting the second component of 
the neutrino mixing matrices in Eqs. (jA-15|) and (jA-16|) : 

w^=E' j U ej V» j = E'Jul : UW) e-*e(v*VW) . (4-14) 

\ / ej \ / 

In this expression, the first matrix element, U e j, is known to be of order unity from neutrino 
oscillation experiments, 5 ^ while, concerning the second element V^j, there is no reliable 
information at present. But if we assume the seesaw mechanism, we must consider (yj )/j 
to have a very small value, as shown in Eq. (jA-6|) . We may get some rough idea of its order of 
magnitude from the neutrinoless double beta decay, which gives the following upper bound 
for a similar quantity: 15 ^ )7 ' ) 

< 0(1(T 6 ). (4-15) 

The primed sum in this case represents a sum that extends over only the light neutrinos 
(rrij < 10 MeV), so that contributions from heavier neutrinos are ignored in comparison with 
the virtual Majorana neutrino momentum. Here, 8 C and 9' c are, respectively, the Cabibbo- 
Kobayashi-Maskawa mixing angles for the left-handed and right-handed d and s quarks. The 
order of magnitude of XW^ in Eq. (j4-13J) seems to be less than 10~ 6 , although the suffix /i in 
Eq. (j4-14|) is replaced by the suffix e in Eq. (j4-15|) . and there are some quantities related to 
the quark sector. Thus the order of magnitude of r) m seems to be much smaller than 10~ 6 , 
because of the additional factor of w efl h in Eq. ()4-13|) . 
*) Our R(x) is equal to 6F T1 (x) of Ref . 1). 



(A)=A 



Ej U e jV e j (COS 9'J COS 9 c ) 
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Recently, Danneberg et al. 4 -* reported the following energy averaged value for e + in the 
direction 9 = n/2: 

(P T1 (x, 9 = tt/2)> = (6.3 ± 7.7 ± 3.4) ■ 1(T 3 . (4-16) 

This experimental result is of the same order of magnitude as the prediction of the standard 
model, 

P rlM ^ /2H _P„ML_| ,4. 17 ) 
where xq is defined in Eq. (|2-8jl . 

4.3. Transverse polarization perpendicular to the decay plane: T(x) 

A non-zero value of T(x) implies the existence of a non-zero Majorana CP violation 
phase in our model. It is expressed as 



T{x) = 2e m ^/x 2 -x 2 r] mi (l-r 2 ) , (4-18) 
where the parameter r] mi is defined as follows, using Eqs. (|B-7|) and (|B-25Jl : 

Vmi = e m {^-^j = £ m (~^~) w^h)- (4-19) 

There is no corresponding term in either the standard model or our model for the Dirac 
neutrino.*-* The parameter i] m i is obtained by taking the imaginary part instead of the real 
part in Eq. (J4-13|) . Therefore, rj m i is proportional to the sin term of the CP violating phases 
appearing in the lepton mixing matrices. As we can imagine from Eq. (j4-15|) . we cannot 
expect to measure i] mi in practice, because its value seems to be too small, as predicted by 
our model (for details, see §3 of Ref. 9)). 

Recently, Danneberg et al. 4 -* reported the following energy averaged value for e + in the 
direction 9 = n/2: 

(P T2 (x, 9 = tt/2)) = (-3.7 ± 7.7 ± 3.4) • 1(T 3 . (4-20) 

A smaller value of (Pt2{%, 9 = tt/2)) is expected if the neutrino is of the Majorana type and 
the CP- violating phase exists. 



*' Our T(x) is equal to 6^7-2(2;) of Ref. 1). Their Ft2(x) includes the CP violating term even in the 
massless Dirac neutrino case, because it conies from the interference between the (V ± A) and (S ± P) ( or 
T ) interactions. 
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§5. Concluding remarks 



It was shown in JJJthat the Michel parameterization, 1 ) which has been used by an ex- 
perimental group, is a special case of the more general form to investigate the deviation 
from the standard model. We propose a new parameterization that directly represents devi- 
ations from the standard model. In general, there is the freedom to choose the normalization 
factor A n £ in Eq. (J3-llj) . Then, we can find the most effective x dependence of the term 
including that parameter, whose value is determined by analyzing experimental data. The 
SU{2)l x SU{2)r x U{1) model has been used to elucidate this feature in the simplified 
expression given in Eq. f)3-30j) . Concerning the spectrum, it has been confirmed that there is 
no significant deviation from the standard model, except that the mean value of pu deter- 
mined experimentally is slightly larger than 0.75, as shown in Eq. (j3-16j) . although we have 
Pm < 0.75 in our model. 

Also, it was shown in £0]that the polarization of e ± can be expressed in terms of the same 
parameters introduced to analyze the spectrum, in contrast to Ref. 1), in which several new 
parameters are introduced. The predictions for the polarization obtained from the standard 
model are also consistent with the recent experimental results. 

It is an important problem to investigate whether the neutrino is of the Dirac or Majorana 
type, as mentioned in §1. In normal muon decay, there are three theoretically possible 
subjects for this purpose within the framework of gauge theory, as we now discuss. 

The first subject is to measure the transverse polarization of e perpendicular to the 
decay plane, namely T(x) in Eq. (j4-18j) . The reason for this is because this polarization 
exists in neither the standard model nor the massive Dirac neutrino case. However, the 
theoretical estimate of T(x) is very small, as explained below Eq. ()4-19p . We cannot expect 
to obtain any useful information from this measurement within the present experimental 
precision. 

The second subject is the transverse polarization in the decay plane, R(x) appearing 
in Eq. (j4-12|) . In this case, the term associated with t] m characteristic of the Majorana 
neutrino has the larger x dependence, but this r] m itself is also too small, as mentioned 
below Eq. f)4- 13j) . It seems difficult to derive any definite conclusion from this measurement. 

The remaining possibility is to take advantage of the different x dependences of the terms 
including the parameters p c and p m in the energy spectrum Eq. (|3-30|) by comparing the \ 2 
values for the Dirac-type neutrino with those for the Majorana-type neutrino, as mentioned 
below Eq. ()3-16|) . This may provide a test to determine the type of neutrino, although it is 
indirect. 

Finally, let us summarize the general features. In the Dirac neutrino case, there is no 



22 



important effect due to the lepton mixing matrices under Condition (B) in Eq. (|2-2()jl . In 
other words, we can use the theoretical expressions obtained by assuming massless neutrinos. 
In addition, we point out that it is useful to choose a different n for the normalization factor 
A nQ in Eq. (3-6) in order to minimize x 2 . Thus, we can find some constraints on the 
coupling constants (A, rj and k) in principle by combining other information from various 
decay processes. We also note that deviations from the standard model become smaller as 
these coupling constants become smaller. 

In the Majorana neutrino case, it is very difficult to find any deviation from the standard 
model under Condition (B) given in Eqs. (|2-21j) and ()2-22j) . This is because all parameters 
include the small components of the lepton mixing matrices. This feature is independent of 
the values of A, 77 and k. 

Appendix A 

Summary of Various Mixing Matrices 



For the purpose of making this paper self-contained, here we summarize the theoretical 
foundations of this work, even though they have been discussed many times in the literature 
already. Many theoretical gauge models beyond the standard model have been proposed to 
analyze normal muon decay. 16 - 1 Among them, let us consider a model that consists of V — A 
and V + A currents. 

The mass term of leptons in the Lagrangian with n left-handed and n right-handed lepton 
doublets is generally defined by 

C M = -T R M E E L - \ (KF, K)) M ^ ( ^ )c \ + H.c, (A-1) 

where £, v L , and v R are, respectively, the weak eigenstates of the charged leptons and 
left-handed and right-handed neutrinos. Explicitly, we write £ R = (e' R , (jf R , •••), £][ = 
( e £>/4'---)> u l = (veL,v»L,---), and v' R T = {v'^v'^ •••)■ Note also that (^l(r)) c = 
C vil(r) T , where C is the charge conjugation operator. Here Me is the n x n mass matrix 
for charged leptons and M. is the 2n x In neutrino mass matrix defined by 

M = ( Mi MS ) , (A.2) 

\ M D M R ) 

where M D , M L and M R are, respectively, the Dirac-type and left-handed and right-handed 
Majorana type nxn mass matrices for neutrinos. Here, the identity v' r MdVl = (^lYM^^u^) 
has been used. 
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Let us first examine the case in which the Majorana-type mass terms exist. Since Ml 
and Mr are symmetric matrices, ^ M. is also symmetric and can be diagonalized by some 
orthogonal matrix in principle to determine the neutrino masses. 18 > However, we use the 
2n x 2n unitary matrix IA V in order to obtain positive values for masses: 7 ^' 19 ) 

U T V MU V = V v . (A-3) 

Here, V u is a diagonal matrix whose 2n elements represent the masses of the Majorana-type 
neutrinos. Therefore, the weak eigenstates of neutrinos are expressed as superpositions of 
the mass eigenstate Majorana neutrinos Nj as follows: 





U V N L = Z N L = - - J ■ (A-4) 



Thus, we have the 2n mass eigenstate Majorana neutrinos, (Ni) T = (A^, N 2 , ■ ■ ■ ,N n ) and 
{Nn) T = (N n+ i, N n+2 , • ■ ■ ,N 2n )- Here, the n x 2n neutrino mixing matrices U v and V u , 
which are expressed as U v = (UP, UP) and V v = {v!) l \ V^), are introduced as 

2n 2n 

vtL = Y^{U v ) t jN jL and Vt R = Y^{V v ) t jN jR (A-5) 

3=1 3=1 

for the case of the n generations. ^ 

In this scenario, the small masses of the left-handed Majorana-type neutrinos (Nj) are 
naturally explained by the seesaw mechanism under the assumption that the right-handed 
Majorana neutrinos (A n ) have large masses. Thus, the elements of both uP and vj 2 ^ are 
of order one, while U {2) and V u {1) are of order m vD /m uR , which seems to be very small. Here 
the quantities m v u and m vR represent the orders of the matrices Mr, and M R , respectively: 

= 0(1), Ui 2) = 0{m uD /m uR ), 

V™ = 0(m uD /m uR ), \/ {2) = 0(1). (A-6) 

From the unitarity condition for U u , the matrices Ujp, ul 2 \ Vj 1 -* and should satisfy the 
relations 

U?>UV» + = 1, (A-7) 

V m v w + V (2) v m = h (A . 8) 

U^V^ T + U^V^ T = 0, etc. (A9) 
Note that UP and vlP themselves are not unitary. 
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We now note that the reason for the smallness of the quantities ue 2 , vg 2 , w eil and w efJL h 
introduced in Eqs. (j2-21|) and ()2-22|) is that they include, respectively, the elements of the 
neutrino mixing matrix products (uPuj> 2)j ) u , (K (1) K (1)t )#, (Uj> 1) V^ )T ) eit and (u!i 1) Vjt 1)T ) lie . 
In fact, the elements like and Vu in them seem to be of order 77Vd/7?vr, as shown in 
Eq. (jXSJ). 

Next, let us consider the other simplified scenario, in which Ml = Mr = 0. Namely, 
there exists only the Dirac-type neutrino mass matrix Mb- Because there is no theoretical 
restriction on Md to be symmetric, it can be diagonalized by two n x n unitary matrices U v 
and V u as 

V}M D U U = D' v . (A-10) 

Here, D' v is a diagonal matrix whose n elements represent the masses of the Dirac-type 
neutrinos. Therefore, the weak eigenstates of neutrinos in this scenario are expressed as 
superpositions of n mass eigenstate Dirac-type neutrinos Nj, as shown in Eq. (|A-5|) . although 
the upper limit of the sum over j is restricted to n. 

In order to avoid the complication of needing to treat the Dirac and Majorana neutrino 
cases separately, we use the following convention in this paper: In the Dirac neutrino case, 

(2) 

only n mass eigenstate neutrinos N\ exist, while An do not exist. Thus we can set U v > = 
and Vy 2 ^ 1 = 0. Therefore, llj) 1 ^ and Vv are both unitary matrices, and they are of order one. 
This situation can be expressed simply as follows: 

£/« = 0(l), U^=0, 

K (1) = o(i), k (2) = o. (A-ii) 

We have considered the simplified case of introducing the Dirac neutrino fields. However, 
there is another possible scenario. In that scenario, one Dirac neutrino field can be expressed 
as a superposition of two Majorana neutrino fields with equal masses. 20 ^ In this case, one 
Majorana neutrino field is chosen from Ai and the other field from An; that is, we have 
rrij+k = rrij for 1 < j(k) < n. 

As a typical example of gauge theory with left- and right-handed weak gauge bosons, 
we consider the SU{2)i x SU(2)ji x U(l) model. 7 )' 17 ) In this framework, the charged weak 
current interaction is written in the bases of mass eigenstates of the charged leptons Eg and 
the neutrino Nj after diagonalizing the charged lepton and neutrino mass matrices: 



2V2 1,3 

- 0f E^7a(l + TsWliVjtiNjRWg + H.c. (A-12) 
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Here, gi and gn are the real gauge coupling constants for the left- and right-handed weak 
gauge bosons, Wl and Wr, respectively, and Ue and Ve are unitary matrices that diagonalize 
the charged lepton mass matrix Me as V e MeUe = De, in analogy to Eq. (jA-10|) . 

The weak gauge bosons Wl and Wr are expressed in terms of the mass eigenstate gauge 
bosons W\ and W% as 

Wl = W\ cos ( + e iip W2 sin (, (A-13) 
W R = -e~ i<p Wi sin ( + W 2 cos C, (A- 14) 



where ( is a Wl — Wr mixing angle and (p is a CP-violating phase. The phase factor e~ lip 
comes from the VEV of the (|, |, 1) Higgs field through the diagonalization of the mass 
matrix of the charged weak gauge bosons. 

In obtaining the effective current- current interaction, the left- and right-handed Maki- 
Nakagawa-Sakata (MNS) lepton mixing matrices U and V are defined as 

U = U%U V = (ului l \U E U^) , (A-15) 
V = e~^V E V u = e-*" (V E V^\ V ] E V^) . (A- 16) 

It should be noted that, although the lepton mixing matrices U and V themselves are unitary 
in the Dirac neutrino case, the 2n x 2n lepton mixing matrix ^ ^ j is a unitary matrix in 
the Majorana neutrino case. 

There exists some freedom in the treatment of the phase factor e~ lip mentioned above. 
In our treatment, this phase factor is absorbed into V of Eq. (|A-16|) in order to make the 
coupling constants rj and K real in the effective Hamiltonian of Eq. (J2-lj) . Contrastingly, 
Herczeg 17 ) includes it in the definitions of k and r\. In our convention, all of the constants 
Gf, A, r] and k are real: 



2 

gii\ A c + tan 2 C 



V2~ 

A= (A ,| 8) 

\g L ) l + A c tan 2 C' V ; 

(g R \ (l-A c )tanC , , 

K ~ V ~ \gj l + A c tan 2 C' 1 } 

where A c = (Mi/M 2 ) 2 , Mi being the mass of W{. Note that we have A ~ k 2 = rf for 
A c < tan 2 C < 1. 
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Appendix B 

Definitions of Various Coefficients 

Our model given by Eqs. (|2-1|) and (|2-3|) differs from the standard model because it takes 
account of the finite neutrino mass and the V + A current. The complete decay formulas 
obtained by keeping both the neutrino masses (rrij and m^) and the V + A parameters (A 
and k = 1]), are given in Ref. 10) 

Terms proportional to the neutrino masses can be omitted in practice, because the neu- 
trino masses are very small in comparison with the energy scale (W): 

1.89- 1(T 7 , 



w) 



(B-l) 



where m u represents rrij and mj,, and is taken to be 10 eV in order to get a rough idea for 
the magnitude. 

Then, we have the following results for various components of the decay rate:*) 
1 



N(x) = ( - 



P{x) 

Q(x) 

S{x) 

R(x) 
T{x) 



,Xj iXj f 



a + (3x — 2x 2 — Xq) + (6 + — a + ) (2x — x 2 — xfy 

+ 12 k + X (1 — x) + 6 E m A d r Xq (1 — x) 

a_(-l + 2x - rl) + (6_ - aJ)(x - r„) 

+12A;_ [l-x] 
a_(3-2x-r 2 ) + (6_ - a_)(2 - x - r 2 ) 
+I2k_{l-x) 
a + (—x + 2 x 2 — Xg)) + (b + — a + ) (x 2 — Xq) 

— 12 k + (—x + x 2 ) — 2 e m A d r x (1 — x) 



X Our 











I 2,. 







' \/ x^ — Xq A <ij (l — r ) . 



(B-2) 

(B-3) 

(B-4) 

(B-5) 
(B-6) 

(B-7) 



Here, we keep all terms with respect to A, k and 77, whereas only their first order terms are 
kept in Ref. 9). The coefficients in these results are defined as follows: 



a± 



(6±A 2 b), 



(B-8) 



(a ± A a), &- 

The relations between the definitions in this paper and in Ref. 9) are as follows: AN(x) = (x/W)N(e), 
AP(x) = -(xp e /WE)P(e), AQ(x) = (x Pe /WE)Q(e), AS(x) = -(x/W)S(e), AR(x) = -(x/W)R(e), 
AT(x) = -(xp e /WE)T(e). 
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k± = (k± c + e m k± m ), (B-9) 

where*-* 

k± c = ^(K 2 c±r l 2 c), k± m = ^K 2 d±r l 2 dy (B-10) 

All of these coefficients are classified into two groups. One group consists of a, b, c, a, b 

and c, which are common to the Dirac and Majorana neutrino cases. The definitions of 

a, b, a and b are given in Eqs. (j2-13|) and (j2-14j) . while c and c are given by 

c = E'^ 2 ( 2i ^ + I^IW. (B-ll) 

jk 

d = J2'l F Jk(' 2F ^ + l^'H^I 2 - (B-12) 

jk 

The other group is only for the Majorana neutrino case: 

d = T,'l F }k(? F ik + o jk ) Re([/;.[/ efe y; j y M ,), (b-is) 

jk 



d = J^'^i 2 (2F jk + G jk ) ReiV^U;^), (B-14) 

jk 

d r = Y^'^kGjkMUljVekV^U^), (B-15) 
jk 

d i = E ' f j1 2g ^ Mu^v;^). (b-16) 
J* 

Under Condition (A) in Eq. (|2-19jl . we have a = 6 and a = b, and accordingly 

a± = 6±. (B-17) 
In this way, the spectrum and polarization given in Eqs. ()B-2|) - (|B-5|) are simplified. 

B.l. Dirac neutrino case 

According to Condition (B) presented in £ 12.21 the masses of the Dirac-type neutrinos are 
conjectured to be so small that all neutrinos are allowed to be emitted in muon decay. The 
finite neutrino masses give rise to slight deviations from unitarity characteristics, through 
the relations given in Eqs. ()2-15|) and (|216|) . We express these small deviations as follows: 

a = 1 - e a {U 2 U 2 ) => 1, b = 1 - e h {U 2 U 2 ) a, (B-18) 

a = 1 - Sa(V 2 V 2 ) ^ 1, b = 1 - e h (V 2 V 2 ) a, (B-19) 

c = 1 - e c {U 2 V 2 ) 1, c= 1 -e £ (V 2 U 2 ) => 1. (B-20) 



*) Here we have introduced the new quantities k±,d r and di in place of x±,h and h 1 used in Ref. 10). 
The relations between them are k 2 x± = 6(k± c + e m k± m ), h = e m d r and h 1 — e m d.- L . 
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Here, for example, e c (U 2 V 2 ) stands for a factor which is written as a product of r 2 k and 
|^ej| 2 |^ufc| 2 , except in the narrow range near x max . The magnitudes of these six factors are 
in general less than 10~ 13 , as seen from Eq. (|2-8jl . The arrows in the above expressions 
represent the limit taken under both Conditions (A) and (B) stated in §21 The values of 
a, a, b, b, c and c can thus all be considered unity in practice.*) 

B.2. Majorana neutrino case 

In the Majorana neutrino case, it is assumed that there exist neutrinos with both small 
and large masses. Heavy neutrinos are forbidden energetically to be emitted in the muon de- 
cay we consider, and the primed sums are therefore only over light neutrinos. The coefficients 
are expressed as follows: 



a => (l-u- e 2 ) (l -u^ 2 ) 

a => v^ 2 v^ 2 , 



d 



w, 



e/i I j 



b 

b 
c 
d 



d r =>• Re(w efM * w efJ , h ), 



a. 

a, 



- 2 (l 



u 



At )i 



I w eflh | , 



(B-21) 
(B-22) 
(B-23) 
(B-24) 
(B-25) 



The quantities ul and here are defined in Eq. (j2-21|) . and and w^Th are defined in 
Eq. ()2-22j1 . The arrows represent the limit taken under both Conditions (A) and (B) given 
in §2j The quantities on the right-hand sides indicate the magnitudes of the coefficients in 
terms of the lepton mixing matrix elements. **' 
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